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I  Introduction 


The  goal  of  this  project  is  the  development  of  techniques  that  will  permit  integrated 
optical  devices  to  be  used  in  inherently  nonplanar  optical  information  processing 
architectures.  As  a  first  demonstration  of  this  method  an  integrated  optical  synthetic 
aperture  radar  (lO  SAR)  processor  is  being  developed.  This  is  a  collaborative  effort 
between  Caltech  and  Professor  A.  Tanguay’s  group  at  the  University  of  Southern 
California.  The  basic  idea  is  an  extension  of  the  space  and  time  integrating  concept 
that  we  have  previously  used  to  demonstrate  bulk  acoustooptic  processors,  including 
a  SAR  system1,2,  to  an  integrated  optics  implementation.  The  bulk  versions  of  these 
systems  make  use  of  all  three  dimensions  of  the  optical  system  as  well  as  time  to  perform 
multidimensional  signal  processing.  Since  integrated  optics  is  a  planar  technology  we  must 
find  ways  to  perform  spatial  integration  in  the  transverse  dimension  of  the  IO  chip  and 
temporal  integration  along  the  direction  of  propagation.  In  addition  we  need  methods  for 
communicating  information  between  the  guided  modes  of  the  IO  structure  and  the  free 
space  modes  in  the  third  dimension. 


In  a  previous  report  (“ Acoustooptic  Synthetic  Aperture  Radar  Processor”,  Grant 
AFOSR-84-0319)  we  have  described  an  architecture  for  implementing  a  SAR  system 
in  integrated  form.  This  architecture  uses  distributed  holographic  coupling  out  of  the 
integrated  plane  to  address  a  CCD  detector  for  temporal  integration.  In  this  report 
we  describe  how  this  distributed  coupling  may  be  implemented  using  gratings  that  are 
holographically  recorded  in  the  plane  of  a  lithium  niobate  waveguide.  These  gratings  can 
be  recorded  by  interfering  a  guided  mode  with  a  beam  incident  on  the  waveguide  from  the 
top.  Thus  we  have  a  method  for  programming  the  operation  of  the  IO  SAR  chip  through 
the  third  dimension.  This  method  of  recording  gratings  for  coupling  guided  modes  to  the 
outside  can  be  naturally  extended  to  implement  more  general  processors.  For  instance,  we 
can  record  photorefractive  gratings  in  waveguides  which  couple  guided  modes  using  light 
incident  from  outside  the  waveguide.  Based  on  this  method  we  have  designed  an  integrated 
optical  architecture  that  is  capable  of  performing  an  arbitary  two  dimensional  linear 
transformation  with  separable  kernel  (or  equivalently  a  triple  matrix  product).  SAR  and 
2-f)  Fourier  transforms  are  two  important  linear  tranformations  that  can  be  implemented 
in  this  manner.  In  section  II  we  describe  the  photorefractive  IO  architecture.  In  the 
following  two  sections  we  discuss  in  detail  the  two  major  components  of  this  architecture. 
We  describe  in  section  III  a  programmable  vertical  outcoupling  mechanism  and  in  section 
IV  an  implementation  of  an  integrated  vector  matrix  multiplier  that  programmable  with 
control  beams  external  to  the  wavoguiding  plane. 
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II  Integrated  Optical  Image  Processing  Architecture 


A  separable  two  dimensional  linear  transformation  takes  a  function  /(£,»?)  into  a 
function  g(£',  rj ')  according  to 

9(Z\v')  =  f  [  (i) 

Jx  Jy 

In  this  architecture  one  of  the  two  input  coordinates  ( rj )  is  represented  by  time.  Specifically 
lines  of  the  2-D  image  are  sequentially  entered  into  the  system.  We  therefore  replace  the 
continuous  variable  77  with  the  integer  variable  n.  The  variables  £  and  £'  are  represented  by 
the  transverse  spatial  coordinates  of  the  chip  at  the  input  and  output  planes  respectively. 
The  output  variable  is  represented  in  the  longitudinal,  or  ^,  direction  of  the  chip.  The 
above  transformation  is  then  written  as: 


g(x',z)=Y2  f{x,n)hi(x,x,)h2{n,z)dx  (2) 

n  ^  * 


In  figure  1  we  show  the  architecture  for  implementing  two  dimensional  separable 
linear  transformations  in  integrated  optics.  The  vector  f(x,n)  is  input  through  a  channel 
waveguide  array.  Light  from  each  channel  is  collimated  by  the  lens  L\.  The  collimated 
beams  are  diffracted  by  a  matrix  h\(x,z')  of  thick  gratings  associating  the  input  channels 
with  an  array  of  output  channels.  Each  of  these  gratings  specifies  a  unique  interconnection 
between  one  channel  at  the  input  and  one  channel  at  the  output.  Crosstalk  between 
channels  is  suppressed  by  the  Bragg  condition  for  volume  holograms,  kx  -  k}  =  KXJ,  where 
ki  is  the  wavevector  associated  with  light  collimated  form  the  ith  input  channel,  k}  is  the 
wavevector  associated  with  light  focused  into  the  jth  output  channel  and  KX]  is  the  spatial 
frequency  of  the  sinusiodal  hologram  connecting  channel  i  to  channel  j.  In  the  case  of 
integrated  interconnections,  all  of  the  beams  of  interest  are  confined  to  a  single  plane. 
In  this  plane  the  two  wavevectors  satilyiiig  the  Bragg  condition  for  a  grating  at  spatial 
frequency  K  are  unique.  This  is  shown  graphically  in  figure  two  for  a  3  x  3  interconnection 
scheme. 

The  lens  £>2  focuses  the  diffracted  input  light  into  the  output  channels.  The  light  in 
these  channels  is  holographically  coupled  out  of  the  plane  onto  a  (’Cl).  In  the  case  shown 
in  figure  1  the  coupling  efficiency  between  the  channels  and  the  CCD  is  modulated  along 
the  length  of  the  waveguides  using  a  l-l)  series  ol  electrodes  addressed  at.  the  sampling 
frequency  to  elcctrooptically  modulate  the  spatial  Irequenries  ol  the  guided  modes  and 
thus  control  the  Bragg  matching  between  these  modes  and  the  radiation  modes  incident 
on  the  CCD.  In  this  way  the  light  detected  at  each  pixel  on  the  (’Cl)  is  modulated  by 


ha (»,*),  where  z  is  a  spatial  coordinate  along  the  length  of  the  waveguides.  h(n,z)  might 
also  be  implemented  by  placing  a  dynamic  spatial  light  modulator  between  the  waveguides 
and  the  CCD.  After  N  pulses  the  charge  distribution  that  accumulates  on  the  CCD  is 
proportional  to 

g(x',z)  =  Ylf  f(x,n)hi(x,x,)h2(ntz)dx  (3) 

Thus  our  architecture  implements  exactly  transformations  of  the  sort  described  in  equation 

2. 

The  gratings  that  are  necessary  for  outcoupling  light  from  the  waveguide  to  the  CCD 
can  be  recorded  holographically  using  the  photorefractive  effect.  This  is  described  in 
the  following  section.  The  gratings  that  are  needed  to  specify  the  linear  transformation 
at  the  front  end  of  the  processor  are  also  recorded  holographically  by  interfering  in  the 
photorefractive  waveguide  beams  incident  from  above.  The  details  of  how  this  is  done  are 
described  in  the  fourth  section. 


Ill  Channel  to  CCD  Coupling 


In  this  section  we  consider  vertical  outcoupling  from  waveguides  for  out  of  plane 
processing.  In  particular,  we  describe  how  vertical  outcoupling  might  be  implemented  in 
the  10  SAR  architecture  described  in  our  previous  report.  The  10  SAR  architecture  is 
completely  analogous  to  the  more  general  linear  transformation  architecture  except  for  the 
fact  that  in  the  10  SAR  case  the  mask  between  the  waveguide  and  the  CCD  need  not 
be  dynamic  on  the  time  scale  of  the  pulse  repetion  frequency.  In  architectures  for  which 
dynamic  masks  are  required,  the  electrooptic  scheme  described  in  the  previous  section 
might  be  employed. 


Synthetic  aperture  radar  is  an  imaging  radar  system  which  uses  the  relative  motion  of 
the  radar  and  target  to  obtain  crossrange  information  from  doppler  shifts  in  the  returning 
radar  pulses.  Range  is  recovered  in  the  10  SAR  processor  by  spatial  integration  of  returned 
pulses  introduced  to  the  device  through  acoustic  transducers.  This  part  of  the  device  is 
described  in  our  earlier  report. 


Azimuth  is  recovered  in  the  IO  SAR  processor  by  a  temporally  integrated  correlation 
of  the  range  focused  signal  against  a  reference  kernel  stored  in  a  two  dimensional  mask 
covering  a  channel  array  to  CCD  interface.  This  mask  modulates  the  field  coupled  from 
the  channels  to  the  CCD  by 


n*.*)  ^ 1  f  co^bll  ))  (*) 

where  62  is  a  constant.  Note  that  we  have  not  included  in  this  equation  the  carriers 
necessary  to  separate  the  full  complex  image  from  the  bias  terms.  The  system  as 
described  here  allows  the  recovery  of  only  the  real  part  of  the  image.  The  carrier  encoding 
mechanisms  are  described  in  detail  in  references  [l|  and  [2|  for  the  bulk  architecture. 
The  same  methods  apply  to  the  integrated  system.  Since  in  the  IO  SAR  processor 
the  light  in  each  channel  must  be  sampled  many  times  along  the  z  direction,  T(x,z) 
must  be  implemented  such  that  the  channels  are  only  slightly  perturbed  at  any  given 
location.  In  this  section  we  describe  how  T(x,/.)  may  be  implemented  using  holographic  and 
photolithographic  techniques  to  modulate  the  amplitude  <>r  perturbations  in  the  refractive 
index  of  the  waveguiding  plane. 


Grating  couplers  are  commonly  used  in  integrated  optics  to  boost  spatial  frequencies 
associated  with  radiation  modes  to  guided  frequencies'*.  A  prototype  geometry  for  these 
devices  is  shown  in  figure  .'5,  I’he  grating  consists  of  a  periodic,  perturbation  to  the 
effective  permitivity  of  the  waveguide.  'I  he  perturbation  may  be  implemented  by  etching 
physical  rulings  as  in  figure  d  or  by  changing  the  actual  index  of  the  waveguiding  film. 
One  method  of  fabricating  gratings  uses  standard  etching  methods  to  record  interference 


patterns  created  by  two  beam  mixing  in  light  sensitive  films.  Grating  couplers  have  been 
made  in  photoresist  by  this  method  at  USC  as  part  of  this  project. 


The  modulation  T(x,y)  could  be  applied  directly  to  a  grating  coupler  by  using 
photolithographic  masks  or  optical  elements  to  modulate  one  or  both  of  the  writing  beams 
during  grating  formation.  Since  the  the  frequency  and  chirp  rate  of  T(x,y)  depend  on 
potentially  variable  aspects  of  the  radar  system  geometry,  the  utility  of  the  10  SAR 
processor  would  be  greatly  extended  if  it  were  possible  to  write  the  grating  coupler  in  real 
time.  In  light  of  this  fact  and  in  light  of  the  long  established  use  of  LiNb03  in  holographic 
applications,  we  consider  photorefractive  holography  directly  in  the  waveguide  array  to  be 
a  most  promising  means  of  vertical  outcoupling  for  the  IO  SAR  and  other  image  processing 
architectures. 


The  photorefractive  effect  is  a  process  by  which  a  spatial  refractive  index  may  be 

written  in  electro-optic  crystals  by  spatially  inhomogeneous  optical  fields9.  The  basic 

mechanism  may  be  visualized  by  considering  figure  4.  Two  coherent  light  beams  with 

wavevectors  k i  and  k2  create  a  fringe  pattern  of  frequency  (ki  —  k2)  in  their  region 

of  intersection.  Charge  carriers  are  excited  from  local  trap  sites  in  the  photorefractive 

crystal.  These  carriers  migrate  to  low  intensity  regions,  where  they  are  trapped  again.  This 

migration  may  be  driven  by  an  external  field,  by  diffusive  relaxation  of  the  inhomogeneous 

carrier  density,  or  by  inhomogeneity  in  the  k  vector  distribution  of  the  excited  carriers. 

In  LiNbC>3  this  last  mechanism,  the  anomalous  bulk  photovoltaic  effect4,  is  dominant 

in  most  applications.  The  migration  of  charge  carriers  from  regions  of  high  intensity  to 

darker  regions  establishes  a  space  charge  field.  The  space  charge  field  in  turn  modulates  the 

refractive  index  of  the  crystal  via  the  electro-optic  effect,  thus  establishing  in  the  crystal 

a  phase  grating  of  period  . 

(*i -*2) 


The  dielectric  permitivity  perturbation  established  in  a  photorefractive  thin  film 
waveguide  by  two  radiation  modes  may  be  modeled  to  a  first  approximation  as 

Ac  =  h(x,y)cos{kz  4-  4>) 


where  k  =  (fci  -  k2)z-  This  perturbation  leads  to  the  exchange  of  energy  between  modes 
at  frequencies  /?„  and  /?„/  where  -  Ak. 


The  coupling  between  mode  u  and 
approximation  by  the  coupled  equations 

d.\, 

dz 

d  A  ,/< 
dz 


mode  i>'  may  be  described  in  the  adiabatic 


J  *  uu’  .4  |  ■' 


J  ^ !/'  u  Au 
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£  and  )i  are  the  normalized  electric  and  magnetic  field  distributions  of  the  modes  and  Au 
is  the  amplitude  of  the  Vth  mode.  Typically  the  coupled  modes  will  be  the  fundamental 
guided  mode  and  the  Bragg  matched  radiation  mode.  In  the  case  of  the  IO  SAR  processor 
we  are  not  concerned  with  the  exact  radiation  mode  into  which  the  guided  light  is  coupled 
since  the  CCD  is  close  enough  to  the  waveguide  array  that  all  free  space  radiation  modes 
will  be  detected.  It  is  important  of  course  that  the  separation  be  large  enough  that  the 
high  dielectric  constant  of  the  silicon  CCD  does  not  perturb  the  waveguide. 


In  the  geometry  we  are  proposing  the  modulation  T(x,z)  is  written  on  the  much  higher 
frequency  carrier  cos (kz  +  <f>).  Since  the  outcoupling  is  driven  at  the  higher  frequency,  the 
distributed  nature  of  the  coupling  may  be  ignored  so  long  as 


T'(x,z) 

T(x,z) 


<  k 


where  T'(x,z)  =  In  this  case  the  only  adjustment  that  need  be  made  to  the 

kernel  T(x,z)  on  account  of  its  longitudinal  implementation  is  that  the  amplitude  of  the 
perturbartion  must  increase  along  the  guide  to  account  for  the  decay  of  the  guided  mode. 


While  a  number  of  experiments  have  been  performed  using  the  photorefractive 
effect  in  LiNb03  indiffused  waveguides  to  write  holograms  transverse  to  the  propagation 
direction5’0,7,  writing  a  grating  along  the  optical  axis  provides  new  challanges.  In  all 
experiments  published  to  date  the  write  beams  have  been  guided.  It  is  not  feasible  to 
use  guided  write  beams  for  our  purposes  because  it  is  difficult  to  form  gratings  along  the 
optical  axis  with  guided  beams  and  because  using  guided  write  beams  in  a  multi-element 
geometry  would  be  damaging  to  the  properties  of  system  outside  of  the  mask  region.  For 
these  reasons  and  for  the  advantages  in  framing  grating  parameters  offered  by  the  extra 
degrees  of  freedom  available  to  radiation  modes,  we  propose  to  use  these  modes  to  write 
the  grating  for  vertical  outcoupling  in  the  IO  SAR  processor. 


Longitudinal  gratings  are  more  difficult,  than  transverse  gratings  to  implement  in 
LiNb03  waveguides  for  materials  reasons.  The  key  mechanism  driving  excited  carrier 
migration  in  LiNbO,{  is  the  bulk  photovoltaic  effect.  The  principal  direction  of  charge 
migration  under  this  effect  is  along  the  c  axis.  For  gratings  written  along  this  axis  the 
applicable  electro-optic,  coefficients  are  r\.>  and  r33.  r:>.>  is  about  three  times  ri3.  If  the 
read  mode  for  the  grating  propagates  transverse  to  /.  the  r33  coefficient,  may  be  read 
by  a  mode  with  the  extraordinary  polarization.  If  the  read  mode  is  longitudinal  then 
the  extraordinary  index  is  not  available  and  the  r  |  -.  coefficient  must  be  used.  This  is 
problematic  both  because  rj3  <  r33  and  because  it.  is  not  desirable  for  the  guided  mode 
to  use  the  ordinary  polarization  for  reasons  to  be  described  below.  One  W’ay  in  which 


these  difficulties  may  be  addressed  is  by  writing  the  grating  along  100  and  using  the  r51 
coefficient,  but  in  this  case  the  grating  must  be  written  by  carrier  diffusion  which  results 
in  low  sensitivity.  Another  possibility  is  to  search  for  crystal  geometries  away  from  the 
crystalographic  axes  that  may  offer  acceptable  compromises. 

Photorefractive  holography  offers  a  flexible  if  challenging  means  of  developing  vertical 
outcoupling  for  two  dimensional  processing.  While  it  is  difficult  to  find  geometries  that 
maximize  both  the  holographic  sensitivity  and  the  waveguiding  characteristics  of  the 
crystal,  it  does  appear  the  geometries  may  be  found  in  which  a  satisfactory  compromise  is 
available. 


IV  Photorefractive  Integrated  Optical  Vector  Matrix  Multiplier 


jj 


In  this  section  we  describe  how  out  of  plane  photorefractive  holography  may  be  used 
to  control  coupling  between  in  plane  guided  modes  and  in  particular  how  photorefractive 
holography  may  be  used  to  implement  an  integrated  vector  matrix  multiplier.  As  discussed 
in  the  previous  section,  photorefractive  gratings  arise  from  the  electro-optical  modulation 
of  the  refractive  index  by  the  space  charge  field  arising  from  the  inhomogeneous  charge 
distribution  induced  by  the  intensity  distribution  of  the  write  beams.  A  particular 
geometry  for  applying  this  effect  in  the  integrated  setting  uses  the  radiation  modes 
Eri{f)  =  Ari£ri{x)e^rl  r  and  Eri{r)  -  A v2£r2{x)e^,-T  to  modulate  the  dielectric 
susceptibility  of  the  waveguiding  substrate  according  to9 

Ax  =  -erE3Ce 

where  e  is  permittivity  tensor  and  r  is  the  linear  electro-optic  tensor.  The  space  charge 
field  is  given  by 

Eac  =  Etf fMiij3\2COs(f}\2  ■  r  +  (j> ) 

Where  M12  is  the  modulation  depth  of  the  intensity  induced  charge  distribution  and 
0i2  —  0i~02-  We  have  assumed  that  M12  1.  The  magnitude  of  Eeff  depends  on  the 

dominant  photocharge  migration  mechanism.  The  three  possible  candiates  are  drift  in  an 
applied  field,  diffusion  and  the  bulk  photovoltaic  effect.  In  the  drift  case  the  space  charge 
field  is  not  linear  in  the  grating  wavevector. 


The  phase  grating  induced  by  the  radiation  modes  couples  the  guided  modes  at  spatial 
frequencies  0gl  and  0g2  such  that  0gi  —  0g2  =  0\2i  where  0i2t  is  the  component  of 
012  which  lies  in  the  plane  of  the  waveguide.  The  component  out  of  the  plane  must 
satisfy  0i2xd  <C  1  where  d  is  the  effective  thickness  of  the  waveguide.  Taking  the 
field  distributions  of  the  co  \pled  modes  to  be  E\{r)  —  Ai{y)£\{x)e](P'r~“t'>  +  c.c.  and 

Ei{f)  =  A2( y)<f2(x)ej('02  +  c.c.  the  coupled  wave  equations  for  A\  and  A2  in  the 

slowly  varying  envelope  approximation  are 


d(l.  lr}  _  r0 

dAn  k 2  f° 

=  j^-EeffMX2(er0l2£)e  /  £i(x)£l(x)e^^'zdx 

dy  2//2  y — oo 

Where  A  =  0\2y  -  ( 0iy  —  02y)-  We  have  assumed  that  the  modulation  depth  is  not 
X  dependent.  If  we  also  assume  that  the  anisotropy  of  the  waveguide  is  small  and  that 
0\2 xd  1  then  for  well  guided  modes 


[  £iix)£i{x)e  l(3vz-Tdx 

J  -  OO 
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where  i/,  is  the  order  of  the  respective  guided  mode.  The  condition  that  the  modes  be 

well  guided  allows  us  to  ignore  the  field  outside  of  the  substrate.  If  tq  =  tq,  the  coupling 

»  2  _* 

efficiency  will  be  large  if  A  <C  ^Eef 

In  the  case  of  vector  matrix  multiplication,  ni«2  gratings  are  used  to  map  nj  input 
frequencies  to  n 2  output  frequencies.  Each  grating  is  associated  with  a  pair  of  radiation 
modes  using  the  rotational  degeneracy  of  the  write  and  read  beams  about  the  grating 
vector.  Consider  the  rotation  a  which  takes  the  input  component  k  1  into  the  radiation 
mode  with  wavevector  K  by  rotation  about  ( ki  -  £2)-  The  equations  aki  =  K  and 
a{k\  —  k 2)  =  (£1  -  £2)  may  be  solved  simultaneously  if  and  only  if  ( K  —  ki)  ■  (Aq  —  k 2)  =  0. 
This  equation  in  concert  with  the  constraints  k^  •  ex  —  0  and  |fc|  =  k0  specify  fc2  to  one 
of  at  most  two  vectors.  This  is  shown  graphically  in  figure  5.  Any  two  vectors  with  end 
points  on  the  arcs  on  the  normal  surface  separated  by  ( ki  — k 2)  which  form  a  plane  normal 
to  the  arcs  may  be  used  to  read  or  to  write  the  grating  with  frequency  (Aq  —  £2)-  However, 
there  need  not  exist  a  single  K  out  of  the  plane  of  &i,fc2>&3  such  that  K  may  be  used 
to  write  both  grating  vectors  shown  in  figure  5.  Thus  it  is  not  possible  to  represent  k  1 
by  a  single  radiation  mode  for  all  712  gratings  associated  with  it.  In  general  as  many  as 
2niri2  points  in  Fourier  space  must  be  available  in  the  radiation  modes  to  write  all  n\ni 
gratings.  One  means  of  presenting  all  these  modes  simultaneously  might  be  to  use  a  second 
volume  hologram  to  store  the  2niri2  beams.  Since  the  intensities  of  all  these  beams  should 
lie  on  the  same  range,  the  modulation  depth  of  each  grating  writen  by  this  approach  will 
be  proportional  to  (nin2)-1.  The  modulation  depths  of  the  individual  gratings  are  also 
proportional  to  (n1n2)_l  for  sequential  writing  if  the  asymmetry  between  the  erase  and 
write  times  is  not  large.  It  is  possible  to  reduce  this  reduction  factor  to  (v/nin2)_I  by 
writing  all  the  grating  with  a  single  reference  beam.  A  large  modulation  depth  is  desirable 
since  the  diffraction  efficiency  of  a  hologram  is  proportional  to  the  square  of  the  modulation 
depth. 


The  architecture  shown  in  figure  6  allows  for  a  simple  writing  geometry  with  a  single 
reference  by  relaxing  the  constraint  of  exact  phase  matching.  A  vector  is  input  through  a 
channel  waveguide  array  integrated  on  a  suitable  substrate.  Light  from  these  channels  is 
collimated  by  integrated  lenses10  before  reading  out  a  matrix  of  holograms  formed  in  a  slab 
waveguiding  region.  The  diffracted  vector  is  refocused  into  an  output  channel  waveguide 
array.  Each  spatial  frequency  in  the  grating  matrix  connecting  the  input  channel  waveguide 
array  with  the  output  array  is  addressed  from  above  the  integrated  plane  by  a  single  pixel 
in  the  input  mask  T.  Light  collimated  from  these  pixels  interferes  with  the  reference  plane 
wave  R  to  create  the  desired  gratings.  The  lenslet  array  LA  may  be  added  to  focus  light 
from  an  extended  beam  onto  the  mask  T  in  order  to  increase  efficiency.  The  modulation 
depth  of  the  (ij)th  grating  is 


Mij 


\J Itj  Ik 
Ik  -t  /V<  /, 


Where  /,y  and  Ik  are  the  intensities  of  the  ( ij)th  write  beam  and  of  the  reference 


respectively.  (/»-,-)  is  the  mean  value  of  the  write  beam  intensities.  The  expected  value 
of  Mij  has  a  maxima  of  ( 2y/(nin2))~ 1  when  (/,,)  =  .  /f?  -  .  The  Fourier  lens  Li  is 

\/(ni ni) 

aligned  so  that  a  pixel  on  its  optical  axis  is  collimated  such  that 

KoOz  =  Kx ;  «oo y  =  Ky  =  0;  /cooz  =  ~KZ 

where  kqo  and  K  are  the  wavevectors  of  the  collimated  beam  and  of  the  reference 
respectively.  The  requirement  that  the  grating  produced  by  the  reference  and  the  light 
from  pixel  ij  form  the  grating  between  the  ith  input  channel  and  the  jth  output  may  be 
expressed  /Cj;  —  K  -  kij.  Since  has  no  component  out  of  the  plane  of  the  waveguide, 
the  ez  component  of  this  equation  is  satisfied  if,  as  we  saw  above, 

(K  —  Kij  ) x  ^ 

In  the  paraxial  approximation  this  requirement  becomes  £  <C  \j ^  where  p  is  the  radius 
of  T.  Having  satisfied  this  constraint,  the  ey  and  ez  components  of  the  phase  matching 
equation  may  be  satisfied  by  properly  placing  the  pixels  in  the  Fourier  plane.  An  example 
of  a  suitable  mask  for  a  four  by  four  device  with  a  400  pm  separation  between  channels  is 
shown  in  figure  7.  We  have  assumed  F=25  cm,  /i  =  /2  =  1  cm.  The  angle  between  the 
writing  beams  is  assumed  to  be  three  degrees. 

We  have  written  photorefractive  gratings  in  single  mode  titanium  indiffused  slab 
waveguides  on  nominally  pure  y  cut  LiNb03.  The  gratings  were  written  using  unguided 
light  at  488  nm  from  an  Ar+  laser.  The  red  HcNe  line  was  used  for  readout.  The  grating 
wave  vectors  were  nearly  parallel  to  the  c  axis  to  make  use  of  the  r33  electrooptic  coefficient 
and  the  photovoltaic  field.  Optical  propagation  was  along  010  for  the  write  beams  and 
along  100  for  the  guided  read  beams.  The  effective  index  for  the  guided  HeNe  mode 
was  2.248,  while  the  extraordinary  index  of  the  substrate  at  488  nm  is  about  2.6.  Thus, 
between  the  write  and  read  beams  there  was  an  expansion  of  the  angle  between  the  beams 
Bragg  matched  to  a  particular  grating  by  a  factor  of  about  1.5.  Small  angles  were  used 
because  our  waveguide  was  not  guiding  along  c.  We  have  obtained  coupling  efficiencies  of 
up  to  4.5%  between  guided  modes  of  interaction  lengths  of  about  1  cm.  We  have  also  been 
sucessful  at  coupling  a  single  guided  mode  into  two  two  diffracted  modes,  though  we  have 
not  yet  implemented  the  mask  scheme  described  above  to  achieve  multibeam  coupling  with 
high  efficiency.  F’hotographs  of  zero  and  first  order  spots  diffracted  from  the  end  of  the 
waveguide  in  the  one  and  two  grating  cases  are  shown  in  figure  8. 


Various  problems  arise  in  the  implementation  of  the  architecture  we  have  proposed. 
Two  specific  problems  concern  the  depletion  of  the  pump  beams  and  the  implementation 
of  the  matrix  in  grating  amplitudes.  Since  the  interaction  region  for  integrated  volume 
holograms  may  be  very  long,  the  coupling  efficiency  of  these  holograms  may  be  quite 
high.  In  the  depleted  pumps  regime  the  accuracy  to  which  the  device  represents  a  true 


vector  matrix  multiplier  will  be  compromised.  The  extent  to  which  this  is  a  problem 
depends  upon  the  application  envisioned.  Since  the  amplitude  of  each  grating  is  linear  in 
modulation  depth,  which  is  a  nonlinear  function  of  the  total  background  intensity,  it  is 
only  possible  in  a  statistical  sense  to  associate  a  given  grating  amplitude  with  a  specific 
transmittance  in  the  writing  mask.  For  small  scale  applications  both  of  these  problems 
may  be  surmounted  by  considering  the  expected  signals  from  each  channel.  For  larger 
scale  applications,  the  application  must  be  tolerant  to  these  nonlinearities.  Optical  neural 
computers11  are  one  such  application,  but  it  seems  unlikely  that  the  device  described  here 
could  be  implemented  on  a  scale  which  would  be  useful  in  neural  computation. 

A  third  problem  arises  from  the  effect  of  titanium  on  the  photorefractive  properties 
of  LiNb03.  Glass  et.  at.  found  that  Ti  doping  increases  the  dark  conductivity  of  LiNbOa, 
thereby  reducing  the  grating  storage  time8.  In  our  waveguides  we  found  grating  lifetimes 
of  several  days  under  continuous  guided  readout  with  about  50  /zW.  One  method  of 
increasing  this  lifetime  might  be  to  use  proton  exchanged  waveguides12.  In  applications 
where  storage  time  is  not  a  concern,  GaAs  or  SBN  might  be  used  as  a  substrate. 
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Figure  S.  Grating  coupler  geometry. 


Figure  6,  Vector  matrix  multiplier  geometry. 


